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SUMMARY 


The information for the determination of the theoretical lift, 
damping -in -roll, and center-of -pressure characteristics for a group 
of plan forms at supersonic speeds has been collected. This information 
is presented in three forms: first, equations; second, figures presenting 

generalized curves in which the quantities plotted are combinations of 
the geometric and aerodynamic parameters; and third, figures presenting 
curves for representative groups of specific configurations showing 
variation of the derivatives with Mach number. From this group several 
observations have been made as to the effect of taper ratio, angle of 
sweep, and aspect ratio on the values of the derivatives. 


INTRODUCTION 


The linearized theory for supersonic flow has been used by various 
investigators to develop expressions for the stability derivatives. 

The information is scattered, however, through a number of reports, not 
all of which use the same notation. A few of the more important 
derivatives are therefore collected herein for future comparison with 
experimental work and to aid designers by presenting the information in 
convenient form. 

In order to make the paper more useful for purposes of comparison 
and design work, calculations were made for related groups of specific 
configurations. The results of these calculations are presented as 
plots of lift-curve slope, damping in roll, and center of pressure 
against Mach number. In addition, all available information is presented 
in the form of generalized figures in which the quantities plotted 
are combinations of the geometric and aerodynamic parameters. 
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SYMBOLS 


7 , 


P 

V 

M 

P 

a 


0 = \]u 2 - 1 


P 

b 

c 


rectangular coordinates (fig. l) 
angular velocity about x-axis (fig. l) 
flight speed 

stream Mach number ( Spee4 sound) 
Mach angle 

angle of attack (fig. l) 

density of air 
span 

local chord 
root chord 
tip chord 


mean aerodynamic chord 


rib/2 


V 


c^dy 


do 


X 

s 

A 

A 
At 

m 
“t 

»^-T 
“ 'AtTTxT 


taper ratio (c-t/c r ) 
wing area 

aspect ratio (b^/s) 

angle of sweep of leading edge 

angle of sweep of trailing edge 

cotangent A 
cotangent A-j. 
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n = = 1 - (1 - X)o) 

“t 

k = \/r 


- pv 

complete elliptic integral of second kind with modulus k 


/ n */2 \ 

yl - l^sin^z dz 

\J0 ) 


F’(Pm) 

E"(M 

I((3m) = 

L 

Z 

M' 


complete elliptic integral of first kind with modulus k 

( M? 


\l yr 


dz 


k^sin^s 


1 

E’(|3m) 


2(1 - P g m2) 


'm 


(2 - p2m2)E'(Pm) - 02m2]r '(pm) 
lift 

rolling moment 
pitching moment about apex 
lift coefficient ( ■=— ^ — 

rolling -moment coefficient 

pitching -moment coefficient 


§pV 2 Sb y 


M 1 


a 


L oa Ja—O 


dC-[ 

da 

dC- 




ipTSSo, 


°”a- 


fSE 

c 


dC 


2k. 

_2V 


►0 


m 


da 


a — >0 


longitudinal distance between leading edge of mean 

aerodynamic chord and center of pressure due to lift, 
fraction of mean aerodynamic chord 
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longitudinal distance between wing apex and center of 


pressure due to lift, fraction of root chord 



Subscripts: 

cp center of pressure 

1,2 regions of wing panel in sketch in appendix 


PRESENTATION OF RESULTS 

Expressions for various stability derivatives have been obtained 
by several investigators (references 1 to 13) with the use of the 
linearized theory for supersonic flow. The equations for dete r mini n g 
the theoretical lift, damping-in-roll, and center-of -pressure character- 
istics of various plan forms have been collected and are given in the 
appendix with a discussion of the range of applicability for each 
equation. All wings are considered to be of zero thickness and the 
equations are given with respect to axes having their origin at the 
leading edge with the x-axis coinciding with the root-chord line. 

Table I lists the various specific wixgs considered according to 
their geometric characteristics. The equations and figure numbers 
given in the table serve as an index to the information presented in 
the paper. References from which the material was compiled are also 
listed in the table. The types of plan forms for which equations are 
given are shown in figure 2. Equations for a few plan forms not 
considered herein may be found in reference 13 . Figures 3( a ) to 3(f) 
illustrate the specific configurations considered. Figure 4 shows 
the variation of the elliptic-integral factors E"(Pm) and l(3m) 
with 3m and is to be used with the equations for the triangular, 
notched triangular, and swept tapered wings in the subsonic -leading - 
edge case ( Pm < l) . 

Computational data obtained from the equations are presented in 
two forms. Generalized curves (figs. 5 to 10), in which the quantities 
plotted are combinations of the geometric and aerodynamic parameters, 
permit a rapid estimation of the value of a derivative for any 
specific wing whose general plan form is considered herein. Curves of 
the lift -curve slope, damping -in -roll, and center-of -pressure 
characteristics plotted against Mach number are presented in 
figures 11 to 1 6 for the specific configurations listed in table I. 
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DISCUSSION 


From the variation of the derivatives (figs. 5 to 1 6 ) as predicted 
by linearized theory, several observations have been made as to the 
effect of sweep angle, taper ratio, and aspect ratio. The following 
effects were noted between Mach numbers of 1 and 2: 

(l) For wings with subsonic leading edges, increase of sweep angle 
decreases the values of and -Cj . At speeds at which the Mach 

line lies behind the leading edge, however, the lift-curve slope and 
the damping-in -roll values are greater with increased sweep. With 
increased angle of sweep, the center of pressure of a notched triangular 
wing moves rearward with respect to both the leading edge of the mean 

5! x 

aerodynamic chord — and the apex of the wing — 

c c r 


(2) Because any change in taper ratio affecth other wing parameters, 
the combined effects of taper ratio and sweep angle were noted with 
aspect ratio remaining constant. For a wing with a supersonic leading 
edge, an increase in taper ratio (ratio of tip chord to root chord) 
decreases the value of the lift-curve slope but increases the amount of 
damping in roll. 


(3) Increase in aspect ratio increases values of the derivatives 
and causes the center of pressure of the wing to move 


An exception is the triangular wing with a supersonic 


Ct and 
-‘-'a 

rearward . 

leading edge where the derivatives and -C^ are functions of 

Mach number alone and the center of pressure is constant at the two- 
thirds -chord point. 


CONCLUDING REMARKS 


A group of 52 specific wings representing various plan forms has 
been considered. The variations of the theoretical lift, damping-in- 
roll, and center -of -pressure characteristics of these wings with Mach 
number are presented. The equations from which these values were 
obtained are also presented. 

From the figures presented, the effect of sweep, taper ratio, and 
aspect ratio may be noted. The lift-curve slope C^ decreased 

CG 

with increase in sweep angle when the leading edge was subsonic and 
increased when it was supersonic. Increased aspect ratio increased 
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the lift-curve slope. In the case of the triangular wing with a 
supersonic leading edge, however, C L was independent of changes in 

sweep and aspect ratio. 


The amount of damping in roll -C 7 was increased with increase of 

P 

taper ratio, sweep angle, and aspect ratio except in the case of wings 
with subsonic leading edges in which an increase in sweep decreased the 
value of -Cjp. Again an exception was noted in the case of the 

triangular wing with a supersonic leading edge for which -Cj was not 


affected by changes in these parameters. The position of the center of 
pressure moved rearward with increase in sweep angle and aspect ratio, 
except in the case of the triangular wing in which it was constant at 
the two -thirds -chord point. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., August 3, 19^9 
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A group of rectangular wings, illustrated in figure 3(a), is 

x c ' 

considered. Variation of PC-^, ~ c > arL ^ with BA is shovn 

r i__ 

in figure 5 and specific curves of Cr , -2£ and -C7 plotted 

■“a. c r ‘'P 

against Mach number are shown in figure 11. Equations (3) to (5) 
are applicable for a range where the Mach line from a tip intersects 

the trailing edge (PA > l) . Equations (l) and (2) for C T and — ^ 

a c r 

are applicable when a Mach line from the tip intersects the opposite 
edge between the midpoint of the side and the intersection of the side 

edge with the trailing edge ^1 > PA j. 

Unswept Tapered Wings 

The equations for and for unswept tapered wings with 

supersonic leading edges which were obtained from reference 5 and an 
unpublished analysis at the Langley Laboratory are as follows : 

Limiting condition: PA - 2 

+ PA)2 (gfc ■ EA)2 . (0 

La ,t3A(P 2 m2 - 1) 2PA(Pm - l) «PA(p2 m 2 - 1)3/2 


_(2Pm + PA)^ P 2 m 2 (pV l4 ' - 4p 2 m 2 + 8)cos’ ] 
48*p3A3(p2 m 2 - i)3 (p2 m 2 - l) 1 / 2 


3Vj; - 10P 2 m2 - 8 
3 


4pA(pm - 1) 


48p2A2(Pm - l)2 


12 pm - PA) (Pm - 2) 

3PA(Pm - 1) + 1 
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The variation of PCt and -PCi with PA for various taper 
ratios is shown in figure 6. The curves of and -Cj , for a group 

of unswept tapered wings (fig. 3("b)), plotted against Mach number are 
shown in figures 12(a) to 12 (d). Equations (6) and (j) are applicable 
when the Mach line from the apex intersects the trailing edge (PA - 2). 


Triangular Wings 

X 

The equations for and for triangular wings with 


subsonic or supersonic leading edges which were obtained from 
references 6, 7, 8, and 9 are as follows: 

(a) Subsonic leading edge 

Limiting condition: PA ^ 4 


pc La . 4f5E»(M 


*CP 


^ = 


2 

3 


-^l(Pm) 


(b) Supersonic leading edge 
Limiting condition: PA = 4 


(8) 

(9) 

( 10 ) 




( 11 ) 


x cp _ 2 


( 12 ) 


PC, 


1 

3 


X 

The variation of PC;L a , and with PA is shown in 


figure 7. The variation of C;L a , -2E, and with Mach number is 


( 13 ) 
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shown in figure 13 for the configurations shown in figure 3(c). 
Equations (8) to (10 ) are applicable when the Mach line from the apex 
lies ahead of the leading edge (pm < l). Equations (ll) to (13) apply 
when the Mach lines from the apex lie behind the leading edge (Pm > l). 


Notched Triangular Wings 


The 


equations for Cq, a > — and for notched triangular wings 


with subsonic or supersonic leading edges follow. They were obtained 
from references 10 and 11. No equation was given for Ci for wings 
with supersonic leading edges. 

(a) Subsonic leading edge 
Limiting conditions: Pm =* 1 

Pm > |n| 

PC L = PAE'^Pm)-^ 1 ' cog~ 1 (-N) + N(1 - N 2 ) 1 / 2 ] 


(1 + K) 


x cp _ N(4 - N 2 ) + (2 + N 2 )(l - Ngr^cos^f-N) 

C r 3(1 - N 2 ) £n + (1 - N 2 )‘ 1 / 2 cos' 1 (-N) 


PC. 


_ -PAl(pm) (1 - N) 


1/2 


48 


(1 + N) 




(14) 


(15) 


3(4N 2 + 1)co3 _1 (-N) + N(2N2 + 13) (1 - N 2 ) 1 / 2 

(16) 


(b) Supersonic leading edge 


Limiting conditions: Pm = 1 


ih^-i 


pc 


La 


8Pm 

0O3 ' 1 (ta) 

jt(N + 1) 

Oft . 2 ->p / 2 


N cos' 


(P 2 m 2 - lW 2 


(il 


(17) 
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(2P 2 m 2 - N 2 - N^cos" 1 ^ j N(3 - ^) CO£f |“ 


^OE (1 - ng)(l 

P 2 m 2 - N 2 ) 3/2 

‘ (1 - N 2 ) ( P 2 m 2 - l)!/2 (p2m 2 - N 2 ) 

G r " 

3 

-cos'll 

(z£) 

[Ml 

1 WcOS~-^T" 

Pm 

( 18 ) 

I(p2m 2 - N 2 )l/2 T (p 2 m 2 _ ±) l/2 


The variation of (3 Ct , and with 3m for various 

•^(X C -p c p 

x 3c 

values of PA is given in figure 8. The parameters Cp^, ~c^* 
and -C i for the group of wings illustrated in figure 3(d), plotted 

Jr 

against Mach number, are shown in figures lh(a), l^b), and lh(c) . 
Equations (lM to (l6) are applicable for a range of Mach numbers for 
which the Mach line from the apex lies ahead of the leading edge (Pm < l) 
and the Mach line from the trailing edge of the root chord lies behind 
the trailing edge (Pm ^ W) . When this latter condition is violated a 
value of the derivative is obtained which has the significance of an 
upper limit (reference 10). This region (Pm < N) is indicated by a 

jx; 

dashed curve. Equations (17) and (l8) for Cp j and — respectively, 

a Cj, y 

are applicable when the Mach line from the apex intersects the trailing 

edge of the wing. The relation between — ^ and — is given by the 

C p c 

expression 

X cp 3 / x cp . A \ 

c 2\c r 12m / 


Swept Untapered Wings 

The equations for Cp, a and for the case in which the Mach 

line from the apex lies on the leading edge (Pm = l) were obtained by 
modifying equations in reference 12 in the following manner: 
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For the lift case. 


°iu - k{°h. + %) 

where subscripts 1 and 2 refer to the regions of the wing panel shown 
in the following sketch: 



If the taper ratio is 1 and the Mach line coincides with the leading 
edge (fim = l), equation (9) of reference 12 becomes 


'Ll 


8m2c r 2a 

rtS 


dv 


(1 - v) 2 \t - V 2 


8b 2 
n S 


a 


ni 


dv 


b-mc 

t/b+mc 


(1 + v)‘ 


•■r 


0>) 


and equation (12) of reference 12 becomes 




2 <y dx dy 
mx + y 


(c) 


Performing the operations indicated in equations (a), (b), and (c) 
yields an expression for in terms of 3 and A. Similarly, 

modifying equation (19) of reference 12 in the preceding manner gives 


°*1 = 


-8n/*c r ^p 


n b-mc 

b+mc 

Jo 


ygdv 

(1 - V)^y/l - v2 


8b3p 

3VS 


v 2 dv 

(1 + v)^\jl - 


3SbV 


(d) 
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and equation (23) of reference 12 "becomes 


'l 2 ~ 3itVSb 



y(3mx - y - b 


h 


~ 7 


/2(mx + y) 


dx dy 


(e) 


The amount of damping contributed by the tip section, however, is seen 
to he negligible and consequently the expression derived for this region 
has been omitted from the expression for C^ . Performing the operations 

indicated in equation (d) and taking the derivative of Cj with respect 

to give the equation for C 7 . 

2 V fc p 


The equations for and derived are as follows: 


Condition of applicability: 


= 1 


PC L = — 
■ Lj a « 


(S5V 3 2 


6 / \ 2 
■ 3® + 84 (* + m) 003 


-1 BA - 1 
PA + 1 


(19) 


3C i = -y 2 
6 P 3 n 




7 15PA 15 p2 A 2 yp3 A 3y 105P 3 A 3 


( 20 ) 


The equations for Cl and C| for a swept untapered wing for the 

u, p 

case in which the Mach line from the apex crosses the trailing 

edge ^ 0 > — + m) were obtained from an unpublished analysis by Mr. Sidney 
M. Harmon of the Langley Stability Research Division and are as follows: 

Limiting condition: 0 > f + — 

A m 




k 

02m 2 (0 2 m 2 - 2)sin -1 ~ 

pm 

- 1 

rt(02m 2 - l) 


££ 


03m3 


jt\/02m2 - 1 2 ( 02m 2 - 1 ) 


+ 0A0m 


( 21 ) 
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_ p5m^(-2p3m 3 + + 2Pm - 9) P^m 3 

Zp 12P 3 A 3 (Pm - l) 2 (P 2 m 2 -1)3/2 3p2 A 2(p2 n 2 _ i)3/2 


p3m3 -2 Pm gV^gVj; - 10P 2 m 2 - 8) 

PA(P 2 m 2 - l) 3 / 2 3^2TI 9«P 3 A 3 ( P 2 m 2 - l) 3 


P^m^ ( p^^ - 4p\A + 8p 2 m 2 )co3 _ -kr- 

+ & ( 22 ) 

3«P 3 A 3 (P 2 m 2 - 1 )^/ 2 


The variation of PC-^ and -PC^ with Pm for various values 

of PA is shown in figure 9- No equation is available for the case in 
which the Mach line from the apex cuts the side edge of the wing; 
consequently the points obtained from equations (19) and (20) have been 
faired to the curves obtained from equations (2l) and (22). These 
regions are indicated by dashed lines. No equations were available 


for 


f2E 


for swept untapered wings. Curves of C^, and 


- 0 \ , 

— u, ‘■p 

with Mach number, for the swept untapered wings shown in figure 3(e) 
are presented in figures 15(a), 15 (b), and 15 (c). 


varying 


Swept Tapered Wings 


Equations for 


Cl and 
L a 



for swept tapered wings obtained 


from reference 12 and from the previously mentioned unpublished analysis 
follow (no equation is given for. Cj^ for the case of a wing with a 

supersonic leading edge) : 
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The variation of 


(3C^ and. -3Ci with 3ta. for various values 


of PA for swept tapered wings of taper ratio 0.5 is shown in figure 10. 
Curves of C L a and -C^ plotted against Mach number are shown in 

figure 16 for the configurations of figure 3(f)* Equations (23) 

(for Ct ) (24) (for -C i_) are applicable for the case in which the 

CO ±r 


Mach line from the apex lies ahead of the leading edge and the Mach line 
from the trailing edge of the root chord lies behind the trailing edge. 
Violation of the latter condition will give values of the derivative 
which have the significance of an upper limit (reference 12) . These 
regions are indicated by the dashed parts of the curves. Equation (25) 
(for Ct ) is applicable when the Mach line from the apex intersects the 

"CL 

trailing edge ^ 3 > ^ ~ n ^ . The two parts of the curves obtained 

from the equations for Ctt have been faired together over the region 
where the Mach line from the tip crosses the side edge of the wing. 
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Rectangu/ar Unwept foperecJ 





(a) Rectangular wings. 

Figure 3»— Specific configurations for which curves are presented. 
Numbers within plan forms correspond to wing numbers in table I. 
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(b) Concluded. 


Figure 3 » — Continued. 
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Figure 11. _ Concluded. 
























(a) Angle of sweep, 30 - 


Figure 1 6 .- Variation of and 

tapered wings. 


_ C •, with Mach number for swept 

P 

Taper ratio, 0.5- 








(c) Angle of sweep, 60°. 
Figure 16.— Concluded. 
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